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Abstract 

A periodically inhomogeneous Schrodinger equation is considered. The inhomo- 
geneity is reflected through a non-uniform coefficient of the linear and non-linear 
term in the equation. Due to the periodic inhomogeneity of the linear term, the 
system may admit spectral bands. When the oscillation frequency of a localized 
solution resides in one of the finite band gaps, the solution is a gap soliton, char- 
acterized by the presence of infinitely many zeros in the spatial profile of the 
soliton. Recently, how to construct such gap solitons through a composite phase 
portrait is shown. By exploiting the phase-space method and combining it with 
the application of a topological argument, it is shown that the instability of a 
gap soliton can be described by the phase portrait of the solution. Surface gap 
solitons at the interface between a periodic inhomogeneous and a homogeneous 
medium are also discussed. Numerical calculations are presented accompanying 
the analytical results. 



1. Introduction 

A homogeneous nonlinear system may admit a localized solutions with a nat- 
ural frequency residing in the first (semi-infinite) band-gap of the corresponding 
linear system. When there is a periodic non-uniformity in the linear system, 
additional finite band-gaps will be formed and the nonlinear system will admit 
a novel type of solitons known as the gap solitons [3]. One main characteristic 
of a gap soliton is the infinitely many number of zeros in the profile of the so- 
lution, inheriting a characteristic of Bloch waves. Gap solitons are intensively 
studied among others in nonlinear optics [2 and Bose-Einstein condensates [18] . 
Several reports on the experimental observation of gap solitons in the fields in 
the one-dimensional setting are, e.g., [6) [TQ| fTT| [23l 124) |3T]. 

Depending on particular underlying assumptions and specific limits, gap soli- 
tons have been studied analytically through several different approaches. The 
first theoretical approach is through the coupled-mode theory, which is based 
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on a decomposition of the wave field into a forward and backward propagating 
wave [33] [6j [8]. The applicability and justification of the method can be seen 
in [12] [26] [27]. The stability of gap solitons in this approach have been stud- 
ied analytically in [22] [7J |4] . The second formal approximation to gap soliton 
is through the so-called tight-binding approximation, which leads to a discrete 
nonlinear Schrodinger equation (DNLS) [17 j. In this approach, a gap soliton can 
be related to the 'ordinary' soliton through the so-called staggering transforma- 
tion. The existence and the stability of discrete solitons in the uncoupled limit 
of this approach has been discussed in [25] . The third analysis of gap solitons is 
based on the approximation when the eigenfrequency of the localized modes is 
close to one of the edges of the finite band gaps [131114 , 28, 34]. In this case, the 
envelope of the gap solitons is described by the nonlinear Schrodinger equation. 
It is shown in [28 that gap solitons at least suffer from an oscillatory instability 
because gap solitons possess internal modes. 

Relatively recently, another analytical method was proposed by Kominis 
et al. [19] [20] [21], employing a phase space method for the construction of an 
analytical solitary wave. Even though the method is rather limited to piecewise- 
constant coefficients, it was shown that the method is effective in obtaining 
various types of localized modes belonging to gap solitons. For that new method, 
the stability result was so far only obtained through numerical simulations. 

The phase-space method proposed in [19] [20] [21] is similar to that used 
in our recent work [16], where it was shown that the profile of a solution in 
the phase-space can be used to describe its instability. The method was based 
on the topological argument developed in [15 . Here, we propose to apply a 
similar method to determine the stability of gap solitons obtained through the 
phase-space method [19 , 20[[2T]. Despite the similarity in the proposed method 
in investigating the instability of gap solitons, the problem is nontrivial. The 
topological argument in [15] is so far immediately applicable to nonlinear sys- 
tems with finite inhomogeneity (see [16] and references therein). By specifically 
constructing the solutions, we show that the argument is also useful to study 
gap solitons. In addition to inhomogeneities occupying the infinite domain, the 
so-called surface gap solitons sitting at the interface between inhomogeneities 
in the semi-infinite domain and a homogeneous region [31] [32] [21] will also be 
studied. Our result will complement the numerical results on the stability of 
surface gap solitons recently studied, e.g., in [9] [5]. 

The paper is outlined as the following. In Section 2, the governing equations 
are discussed and the corresponding linear eigenvalue problem is derived. The 
construction of gap solitons using the phase-space method is briefly explained. 
The instability of gap solitons is studied analytically in Section 3 using the 
topological argument. In Section 4, the linear eigenvalue problem for several 
gap solitons is solved numerically, where an agreement between the analytical 
results presented in the previous section is obtained. In the same section, the 
instability of surface gap solitons is also discussed. We conclude the paper in 
Section 5. 
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2. Mathematical model 

We consider the following governing system of differential equations 
i*t + *xx + l^l 2 ^ = VV x eU := R \ £/> 



where the 'outer' equation has focusing type nonlinearity, the 'inner' equation 
can be defocusing (77 > 0) or linear (77 = 0), and Uo,Ui are disjoint sets of 
intervals to be specified later. 

To study standing waves of 0, we pass to a rotating frame and consider 
solutions of the form ty(x,t) = e~ lUJt ip(x^t). We then have 



#t + ipxx + M 2 ^ = (V -uj)ijj x e U 
#t + ^ - r/l^l 2 ^ = x £ XJj. 



2.,, _rr. (2) 



Standing wave solutions of ([!]) will be steady-state solutions to ([2|. We consider 
real, t independent solutions u{x) to the ODE: 



(3) 



To obtain solutions that decay to as x — > ±00, the condition that V — > 
is required, with G M+. We will also require that u x —> as x — > ±00. To 
establish the instability of a standing wave solution we linearize Q about a 
solution to ([3]). Writing ip = u(x) + e + is{x))e xt + (r(x)* + zs(:r)*)e A *) 

and retaining terms linear in e leads to the eigenvalue problem 

r \ ( D_ \ ( r \ M / r \ (4) 



s i y — D + J \ s J \ s 
where the linear operators D+ and D_ are defined as 

D+= JJ-(^-^) + 3^ 2 , xeUo, (5) 
^2 +cj - 377W 2 , x eUi, 

D = j^-(V-uj) + u\ xeU , (6) 
J^+u-rju 2 , xeUi. 

It is then clear that the presence of an eigenvalue of M with positive real part 
implies instability. 

In [19] a gap soliton was constructed via a method of superimposing the 
phase portraits of the 'outer' system: 

u x = y, y x = (V -u)u-u 3 , (7) 

and the 'inner' one: 

u x = y, y x = -um + r]u 3 . (8) 
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We can view the composite picture as a single, non- autonomous system with 
phase plane given by: 

u x y, 

_ f (V -uj)u-u 3 , x e U , (9) 

~ \ -UJU + T]U 3 , X e Uj. 

In the phase plane of Q, the outer system admits a soliton solution, given 
by the equation: 

y 2 = (V-c)u*-^, (10) 

while solution curves of the inner system are given by 

y^-^ + Tl+C. (11) 

The inner system (|8| admits a heteroclinic orbit in the phase plane given by 
C = cu 2 /2. The solutions we are interested in will travel in the phase plane along 



the homoclinic orbit of the outer system described by (10) and then 'flip' to the 
inner system as x passes through Z7/, and then 'flip' back to the outer system 
along the homoclinic orbit, repeating the process for each of the components of 
XJi (see [19]). 

Let Us be the collection of intervals Us = [0, xq) U (xi, X2) U (#3, x±) — In 
the case of a gap soliton, Uj = -C/5 U C/5, and we have that the number of 
components of Uj is infinite and the x[s are chosen so that the soliton travels 
from (^0,2/0) along the inner system to (— uo,—yo ). This is a key ingredient 
in the construction of the soliton, and will play a large role in establishing 
instability. In [19 , the inner system is linear, and the length of the interval 
(x2k,X2k-i) can be determined as nj^/uj. Here, we do not require that the 
inner system be linear, however we do require that the x[s be chosen so that if 
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i > 1, the soliton travels from (uo,yo) on the homoclinic orbit along the inner 
system to (— uq, —yo), which is also on the homoclinic orbit. 

In Figure [T] we plot an example of a gap soliton of the governing equation 
([I]) for parameter values that will be explained in Section [4j One can notice 
the main characteristic of gap solitons in the plot, which is the infinitely many 
zeros in the soliton profile. 



3. Instability Results 

To show instability of the standing waves, we will show that the matrix M 
from above has a real positive eigenvalue. This is done by applying the main 
theorem of [15]. In [16], systems like ([!]) were considered with Ui = (— 
for some real number L. One can show that the following quantities are well 
defined (see for example [15], and the references therein): 

P = the number of positive eigenvalues of D+ 
Q = the number of positive eigenvalues of D_. 

We then have the following: 

Theorem 1 (|15j). If P — Q ■=£ 0,1, there is a real positive eigenvalue of the 
operator M . 

From Sturm-Liouville theory, P and Q can be determined by considering solu- 
tions of D+v = and D_v = 0, respectively. In fact, they are the number of 
zeros of the associated solution v. Notice that D_v = is actually satisfied by 
the standing wave itself, and that D+v = is the equation of variations of the 
standing wave equation. It follows that: 

Q = the number of zeros of the standing wave u. 

P = the number of zeros of a solution to the variational equation along u. 

(12) 

For gap solitons, it is not immediately clear how to apply Theorem [l] above 
as in this case, both, P and Q — >• oo. The idea presented in this paper is to 
build an approximation to a gap soliton using more and more intervals of Ui for 
which the quantity P — Q remains constant. To this end define 5o = [0, xq) and 
S n = [0,x ) U (x u x 2 ) U (x 3 ,x 4 ) U . . . f> 4n _i,£ 4 J, where (x i: x i+1 ) C U s . Thus 
S n adds two more components for each n. Then we can define U n = —S n U S n , 
and we let f n be a solution to the ODE 



fxx = (V-uj)f-f, xeR\U n , 

fxx = + X G U n . 

Thus for example fo would be the solution to 

fxx = (V-Uj)f-f 3 , \x\>X , 

fxx = -Uf + T]f 3 , \x\ < X , 



(13) 



(14) 
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while fi would be a solution to 

fxx = (V - u)f - / 3 , x £ (-x 4 , -x 3 ) U (-x 2 , -xi) U (-x ,x ) U {x u x 2 ) U (x 3: x A ) 
fxx = -wf + vf 3 , x G (-x 4 , -x 3 ) U (-x 2 , -xi) U (-x , x ) U (xi, x 2 ) U (x 3 , £4). 

(15) 

A gap soliton then can be realized as the limit of successive / n 's (in a variety 
of norms, but in particular in the L 2 and H 1 norms). In Figured we present a 
plot of / n , n = 0, 1, 2, approximating the gap soliton in Figure [l] 
We have the following theorem 

Theorem 2. The quantity P — Q is the same for all fi described above. Thus 
if fo is unstable then so is f n for all n. Further, if /o is unstable, then so is f , 
the gap soliton, corresponding to the limit. 



The key idea is to use the interpretation of P and Q given in (12) as the 
number of zeros of the solution / and the number of zeros of the solution to 
the variational equation along /, for the partial solution defined on (2^,2^+4), 
to the ODE below: 



fxx = (V-CJ)/-/ 3 , X e (Xi+i,Xi+2) U (x i+3 ,X i+4 ) 

fxx = -uf + nf 3 , xG(x,,z i+ i)U(% 2 ,%3). 



The number Q is straight forward to calculate. We make the geometric 
observation as in [16] that P, the number of zeros of a solution to the equation 
of variations along /, can be found by determining the number of times that a 
vector must pass through the vertical as the base point ranges over the entire 



orbit. It turns out that for the solution of (16) defined above, the rotation of a 
vector by the equation of variations is the same (mod 2tt) as if the base point 
had traveled along only the outer homoclinic orbit. 

Example 1. To better illustrate this last point, we first consider the case when 
both the inner system and the outer systems are linear. That is we have the 
following systems of linear, constant coefficient equations 

u\ _ ( 1\ (u 
VJ„ " \(V-«>) Oj [y 
1\ (u 

-LJ J L 



, when x e (x i+1: x i+2 ) U (x i+3: x i+4 ) (17) 
, when x e (xi,x i+ i) U (x i+2 ,x i+3 ) (18) 



The solution to the above equation can be written explicitly. Further, because 
we are in the linear case, we have that the equation of variations along a solution 
is the same as the equation itself ( T7fll8 ). 



Being led by the geometry of the phase plane, we let 3>i(a, b) denote a 
fundamental solution matrix to the equation of variations of the outer system 



of equations (17) along a solution to (17) which travels from point a to point b 



in the phase plane. That is let (u(x),y(x)) be a solution to (17), considered on 
the interval (xj,Xk)- Then set a := (u(xj),y(xj)) and b := (u(xk),y(xk)), and 
define 3>i(a, b) to be a fundamental solution matrix of the equation of variations 
to the outer system, along the path (u(x),y(x)) with x G (xj,Xk)- 
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Figure 3: A sketch of a phase portrait of the partial solution to equation \16\ . The points ai 
correspond to the points (/(a^_i), f x (xi—i)) in the phase plane. 



Similarly, let $2 (a, 6) be a fundamental solution matrix to the equation of 
variations of the inner system (18), along a solution to (18) evolving from point 



a to point b. We denote by ao, ai, &3, the points in the phase plane of (17 18) 



we stop evolving (see Figure 



3 ) , and we let 



be a pair of initial conditions 



where the solutions switches between the two systems, and a 4 the point where 

in the tangent plane to R 2 at the point ao- We have that a solution to the 
equation of variations along the orbit from ao to a\ to a2 to a 3 to a 4 can be 
described as 



$i(a 3 , a 4 )$ 2 (a 2 , a 3 )$i(ai, a 2 )$ 2 (a , a i) 



Co 



It turns out that modulo 27r, 



$1(^3, ^4)^2(^2, a 3 )$i(ai, a 2 )$ 2 (ao, «i) 



= $i(a ,a 4 ) 



(19) 



The equality in equation ( 19 ) can be verified by solving the appropriate systems. 



Another way to see the effect is to consider the following. As the base point 
evolves under equation (T7|[l8) from to a^+i, we can consider the aggregate 



effect of a 3> 7 (a^, a^+i) on a tangent vector ( t ] , as a linear map from 
by simply determining where a tangent vector to a$ gets sent to, when the base 



Co 
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point is at a^+i. That is we are considering $j(a^,a^ + i) as a map between the 
tangent plane of R 2 at the point a$ to the tangent plane of R 2 at the point 
Oj + i. This will give us the total rotation of a tangent vector modulo 2tt as 
we travel from point to point a^+i along the orbit. The key observation is 
to realize that for $2(^0? an d $2(^2? ^3)? this will be negative the identity 
—Id. That is, viewing $2(07, aj+i) j = 0, 2 as a map between tangent spaces 
of R 2 , $ 2 (%,%+i) : T a .R T a . +1 R j = 0,2, we have $ 2 (a j ,a j+1 ) = -Id. 
Moreover, by considering $2(07, ay+i) in this way, we are just measuring the 
effect of rotation by $2(07"? o?"+i) on an initial tangent vector modulo 27r, and 
we have that 

$i(a 3 , a 4 )$ 2 (a 2 , a 3 )$i(ai, a 2 )$ 2 (>o, ai) ^ 
= (-Id) 2 ^ 1 (a 3 ,a 4 )^i(a 1 ,a 2 ) f^j (20) 
= *i(a ,a 4 ) ^ 

where the last equality follows from the facts that clq = — a\, a 2 = — a 3 , the 



outer system of equations (17) is symmetric about the origin, and the group 
property of variational flows. 

We are now ready to state the main lemma used in the proof of theorem [2] 

Lemma 1. Redefine 3>i(a, b) and $2(^,6) as in the above example, but instead 
of using the linear ODE, let them be the fundamental solution matrices to the 
equations of variations along solutions to the inner and outer systems given in 
the nonlinear equation (it 



fxx = (V-Uj)f-f 3 , X e (Xi+i,Xi+ 2 ) U O i+3 ,£ i+4 ), 
fxx = -^/ + ^/ 3 , xe(x i ,Xi+ 1 )U(Xi+2,X i + 3 ). 

Likewise, let aj be defined analogously for the points in the phase plane of the 
nonlinear equation where the orbit switches between the inner and outer systems. 

Also, let ( ° ] be an initial condition to the equation of variations along a 



solution to (16) in the tangent plane to R at ao. Then we have the following: 



$i(a 3 ,a4)$2(a2,a 3 )$ 1 (a 1 ,a 2 )$2(ao,a 1 ) ( £° ) = $i(a ,a 4 ) ( £° ) • (21) 



Proof. The exact same reasoning can be used to prove Lemma[l](the nonlinear 
case), as was used in the example (the linear case). The only difference is that in 
order to determine the aggregate effect of the inner system on an initial tangent 
vector some more care must be taken with the matrices 3>2(&i> a^+i). Write the 
equation of variations to the outer system as 

' ^ - ' ,, n J ( i J ' when x G (x i+ i,x i+2 )U(x i+3 ,x i+4 ), (22) 



v £ J \-3u( + V - uj Oi ^ 
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where u\(x) is the equation satisfying the outer system with linx^ioo u\(x) 
lim^^ioo u[(x) — 0. Write the equation of variations of the inner system as 



D (3^ -co 0) U 



when x G (0^,0^+1) U (0^+2,^+3), (23) 



where i^ 2 , satisfies the appropriate conditions for the orbit. Now here is where 
the appropriate choices of the a^'s must come into play. In the linear case, the 
XiS were chosen so that the length of an interval in Ui was ^!=. Here we choose 
the x^s in Ui so that the length of an interval is such that we will return not 
only to the homoclinic orbit, but also if we leave the homoclinic orbit at the 
point (uq, 2/0)5 we will return to the homoclinic orbit at the point (— uo,— JJo)- 
This allows us to determine the effect of the rotation (modulo 2tt) by the flow 
associated to the equation of variations along the partial orbit (1^2 (x), 1/2 0*0 )• 
In fact, we claim that the exact same is true as in the linear case. If B is the 
linear map from the tangent space at ao and at d2 to the tangent spaces at ai, as 
respectively, then B = —Id. To see this we will write out B in a suitable basis 
vi,V2 of the tangent space at ao- One obvious choice of a basis vector is the 



tangent vector to the inner system. However given equation (23), and the fact 
that along an orbit (uo,yo) — >• (— uo,— yo), this means that if v\ is the vector 
tangent to the inner orbit at ao (or 02), then under B v\ —> —v\. This means 
that B has the form: 

B=(-' (24) 



62,2^ 

where bij are the coefficients of the linear combination of v{ and a suitably 
chosen v^. Now we appeal to two facts about the matrix B which are evident 
from it's definition. The first is that B must be orientation preserving. This is 
an elementary consequence due of the fact that it is the matrix of a flow (see 
for example [29]). This means that 62,2 must be negative. The second fact is 
that since B corresponds to the matrix of the equation of variations traveling 
half way along the periodic orbit given by (1*2 2/2(2?) ) (because we chose our 
x^s so it would be that way), we must have that B 2 = Id. But this means that 
&i,2 = an d ^2,2 = — 1 and the matrix B itself B = —Id. Now we simply repeat 
the computation done in equation (20) and the proof of Lemma [l] is complete. 



We are now ready to complete the proof of theorem [2] . 

Proof of Theorem El Recall that f n as constructed is the solution to the 
ODE ( [l3]) . We let P n and Q n denote the count for f n of P and Q respectively. 
Lemma[l| shows that P n _i = P n + 2 and it is clear that Q n -\ = Q n , and so the 
quantity P n — Q n is the same for all f n , and in particular is equal to P — Q for 
fo. This completes the first part of the proof of theorem [2] 

In order to determine the instability of the limit soliton we must proceed 
topologically using the methods developed in the proof of the main theorem of 

eng. 
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We have already discussed that in H 1 , f n —> f a solution to 
f xx = (V-co)f-f x€R\U!, 



I XX 



Ljf + T]f 3 x ^ C^/- 



Following [15] we can associate to each solution f n a curve j n (x), and to / 
a curve j(x) in A(2) the space of Lagrangian planes in R 4 . 

This is done as follows. Let $2 (x) denote the evolution operator of the 
ODE corresponding to the equation of variations of the ODE (13) along the 
solution f n . Likewise, let <&l+(%) denote the evolution operator of the ODE 
corresponding to the equation of variation of the ODE ( 25 ) along the solution 

/ = lim f n . Thus if ( ] is a pair of initial conditions at x = 0, then for any 

n— )-oo \WQ J 

x G R we have that the evolution of ( V ° ] under the equation of variations along 

V^o J 

/, /„ respectively will be given by (^^ + ^) . ^° ) ' res P ectivel y (tl + (x) • w 

We remark that the initial conditions ( V ° J will be the same for each f n as 
well as for /. 

Again appealing to [15], we can explicitly write the curves j n (%) an d j(x) 
in the space of Lagrangian planes A(2) ~ U(2)/0(n). This is given by 

where 

(x) • w f ' (x) 

e ltn = 2arctan(-^- ) and, # 2 , n = -2arctan(^fA (27) 

( o A)J ' (28) 

01 = 2arctan( t L+ ^ ! '^° ) and, 2 = -2 arctan(4^v)- (29) 

Now we observe that the curves j n (x) and j(x) actually lie on a torus contained 
in A(2). 

It was established in [15 that because f n and / are solutions corresponding 
to homoclinic orbits in the phase plane of equations (13), and (25), the curves 
7 n (x), and j(x) have well defined end points. Let /i-, n , M+,n be the endpoints 
in A(2) of j n (x). That is, let 

lim j n (x) /i_, n and, lim x ^ QO 'y n (x) [j,+ n , 



and 



where 
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and set 

lim j(x) — /i_ and, lim x ^ 00 j(x) = 
x — y — oo 

Further because f n — » f and lemma[TJ we have that /i_ ?n = //_ , and /i+, n = 
for all n. In the previously introduced coordinates on the torus in A(2) this 
means that the limits of d\, n , #2,n are equal to the limits of 0\(x) and 02(x) as 
x — » ±00. Moreover, it is easy to calculate explicitly that 

9\{x) -» 2arctan(v / V - cj) := 0_ and, 6> 2 (x) — 0_ 

as x — » —00. 



0.1 



0.05 
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□ •□ □ ffl is 
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Re(X + ) 



Figure 4: The positive eigenvalues of the operator D+, i.e. A+. One symbol corresponds to 
two different, but very close eigenvalues. 

Still following the outline laid out in [15], we denote by ~ the lift of the point 
(or curve) in the torus embedded in A(2) to its corresponding point in the 
universal cover of the torus, R 2 . We will parametrize the universal covering of 
the torus in the obvious way. Without loss of generality, all of the /i_, n 's and 
/i_ can be lifted to the same point jl- = (6L, —9_). It was shown in [15] that 
for each n, /i +?n lifts to the point /i +5n = (±6L,#_ + (P — Q)2n. Thus lemma 
[l] implies that each /i +5n lifts to the same point /i +) o = (±0_,0_ + 27rfc), 

Next we observe that as f n — ^ / pointwise, j n (x) —> j(x) in the torus inside 
A(2) pointwise, and the compactness of the torus and of A(2), means that the 
end point /i + must lift to the same point in the cover as /i+,o- Thus we have 
that /i + = (±6L,6L +2tt/c). 

Finally, it was shown in [15], that if \k\ 7^ 0, 1, then the corresponding soliton 
underlying the curve 7 is unstable. This completes the proof of theorem [2] 
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Remark 1. The proof of theorem 2 may also be couched in the language of 
fixed end point homotopy classes. There are several ways to define such classes, 
see for example [30] or pQ, and the references therein. In this context theorem [2] 
establishes that the fixed end-point homotopy class of the curve 7 is the same as 
those for j n (x). An immediate consequence of this observation is that in A(2), 
it is possible to deform the curves 7, and 7 n all to the curve 70, in a continuous 
way. 

Remark 2. One can also consider so-called surface gap solitons, and obtain 
exactly the same results as for theorem [2] Mathematically, a surface gap soliton 
is the evolution of the solution to equation ( [25] ) but with the chosen intervals C/j 
replaced by Us, defined earlier. In this case, we consider a sequence of functions 



f n which are solutions to equation (13), but with U n replaced by S n . Then the 



functions f n — » /, a solution to ( |25| ) with the appropriate replacements. Lemma 
[l] holds, as well as theorem [2j and the techniques used in each will be identical. 
Thus if we start with an unstable solution, then the surface gap soliton that we 
obtain in the limit will also be unstable. (See below for a further discussion of 
surface gap solitons). 



4. Numerical solutions and Discussion 

We have solved the time independent equation ([3| numerically, where we 
have used a spectral difference method to approximate the Laplacian u xx . Once 
a solution is obtained, the corresponding eigenvalue problem Q is solved using 
a MATLAB routine. The time dependent equation is integrated numerically 
using a fourth-order Runge-Kutta method. Throughout the paper, we consider 
the parameter values 

V = 1, oj = 0.5. 
First, we study Equation with 

f 1, x e (-ffo,a ), (30) 

\ 0, X E {x 2n+1 ,X 2n+2 ), (-^2n+2,-^2n+l), 

where x = 2, x 2n+1 - x 2n 1, ^2n+2 - x 2n+1 tt/^/u and n = 0, 1, 2, A 

gap soliton for the above periodic inhomogeneity is depicted in Figure [I] 

Theorem [2] implies that to determine the instability of the gap soliton, it 
suffices to determine the instability of the corresponding solution fo shown in 
panel (a,b) of Figure [2] As discussed in [16 , the positive solution f is unstable, 
with P = 2 and Q = 0. We plot A + , i.e. the eigenvalues of the operator D + , 
in Figure [4j As shown in the figure, for fo there are two positive eigenvalues 
of D + , i.e. P = 2. The matrix M in Q for the solution has one pair of real 
eigenvalues [16] in agreement with Theorem [l] 

According to Lemma [l] f n must have the same value of P — Q as /o . In 
the same figure, we obtain that f\ and f 2 respectively has P = 6 and P = 10. 
Considering the fact from Figure [2] that f\ and f 2 respectively has Q = 4 and 
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Q — 8, we indeed obtain that P — Q = 2 for both fi and fa. Using the lemma, 
one will obtain that P — Q = 2 for limn^oo f n . Using Theorem [2j one can 
conclude that the gap soliton in Figure [I] will be unstable. We depict in Figure 
|5|a) the eigenvalue structure of the gap soliton in the complex plane. When the 
corresponding fa of the gap soliton has one pair of real eigenvalues [16] , the gap 
soliton has several pairs of unstable eigenvalues. Nonetheless, one can easily 
notice that there is only one pair of real eigenvalues, similarly to fa [16 j. The 
time dynamics of the solution is shown in panel (b) of the same figure, where a 
typical instability is in the form of the dissociation of the solution. 
Next, we study Equation ([I]) with 



x e (-xo,x ), , 31 x 

X e 2n +l,^2n+2), 



for the same values of x ni n = 0, 1, 2, . . . , as above. The only difference with r] 



defined in Equation (30) is that the present periodic inhomogeneity only occu- 
pies the x > 0-region. In this case, we will have surface gap solitons sitting at 
the interface between a homogeneous and a periodically inhomogeneous region. 
A corresponding surface gap soliton of that in Figure [I] and one of its successive 
approximations fa are shown in Figure [6] The fa approximation of the soliton 
is nothing else but that shown in Figure|2|a). 

Using Theorem[2] and Remark [2j one can expect that in this case P — Q = 2. 
Plotted in Figure [7[a) is the positive eigenvalues of D+, i.e. A + . The positive 
eigenvalue A + of fo is the same as before, which is P = 2. For f\ and fi, from 
Figure [TJa) one can deduce that P = 4 and P = 6, respectively, with Q = 2 
and Q = 4. Hence, the limiting quantity P — Q of the surface gap soliton is 
the same as that of the gap soliton in Figure [l] i.e. P — Q = 2. As expected, 
shown in Figure [TJb) is the eigenvalue structure of the gap soliton, where one 
also obtains one pair of real eigenvalues similarly to the stability the gap soliton 
depicted in Figure [5| a). We plot the time dynamics of the surface gap soliton 
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Figure 7: (a) The positive eigenvalues of D+ for the approximations f\ and $2- Note that 
different from the plot in Figure [J] here each symbol corresponds to one eigenvalue, (b) The 
eigenvalue structure of the surface gap soliton in Figure |6ja). 
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Figure 8: A time dynamics of the surface gap soliton in Figure [6] Shown is the top view of 
\ip(x,t)\ in the (£, x)-plane. 

in Figure [8] 

5. Conclusion 

We have considered a nonlinear Schrodinger equation with periodic inho- 
mogeneity, both in the infinite and semi-infinite domain. Specifically we have 
studied the instability of gap solitons admitted by the system. We have estab- 
lished a proof that if the periodic inhomogeneity is arranged in a particular way, 
such that parts of the solutions belonging to closed trajectories in the phase- 
space have length half the period of the trajectories, then the solitons inherits 
the instability of the corresponding solution with finite inhomogeneity. The an- 
alytical study is based on the application of a topological argument developed 

in na. 

It is natural to extend the study to the case when the solutions are localized, 
but do not tend to the uniform zero solution (see, e.g., [20). The (in)stability 
of such solitons is proposed to be studied in the future using analytical methods 
similar to that presented herein. 
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